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The S shape of the canonical heat capacity is known as a signature of the pairing transition and 
along an isotopic chain it is significantly more pronounced for nuclei with an even number of neutrons 
than with an odd number. Although the heat capacities extracted from experimental level densities 
in 93_98 Mo exhibit a clear S shape, they do not show such an odd-even staggering. To understand the 
underlying physics, we analyze thermal quantities evaluated from the partition function calculated 
using the static-path plus random-phase approximation (SPA+RPA) in a monopole pairing model 
with number-parity projection. The calculated level densities reproduce very well the experimental 
data and they also agree with estimates using the back-shifted Fermi-gas model. We clarify the 
reason why the heat capacities for Mo isotopes do not show odd-even staggering of the S shape. 
We also discuss thermal odd-even mass differences in 94_97 Mo using the three-, four-, and five-point 
formula. These thermal mass differences are regarded as indicators of pairing correlations at finite 
temperature. 

PACS numbers: 21.60.Jz, 21.10.Ma, 05.30.-d 



Level density is one of the basic ingredients required 
for theoretical studies of nuclear structure and reactions. 
Recently, level densities for 93-98 Mo [lj have been ex- 
tracted by the group at the Oslo Cyclotron Laboratory. 
In their paper, the S shape of the canonical heat ca- 
pacity is interpreted as consistent with a pairing phase 
transition with a critical temperature for the quenching 
of pairing correlations of T c ~ 0.7-1.0 MeV. 

Pairing correlations are of special importance in nu- 
clearphysics. The Bardeen-Cooper-Schrieffer (BCS) the- 
ory [J has successfully described the sharp phase transi- 
tion connected to the breakdown of pairing correlations 
for infinite Fermi system of electrons in low-temperature 
superconductors. This sharp phase transition leads to 
a discontinuity of the heat capacity at the critical tem- 
perature, which indicates a second-order phase transi- 
tion. For finite Fermi systems such as the atomic nucleus, 
however, nucleon number fluctuations and thermal and 
quantal fluctuations beyond the mean field become large. 
These fluctuations wash out the discontinuity of the heat 
capacity in the mean-field approximation and give rise 
to an S shape |3( ■ In recent theoretical approaches such 
fluctuations have been taken into account. Examples 
are: the static-path plus random-phase approximation 
(SPA+RPA) [land shell-model Monte-Carlo (SMMC) 
calculations 

It has recently been reported @ that the canonical 
heat capacities extracted from experimental level den- 
sities exhibit S shapes with a peak around the critical 
temperature. The S shape is well correlated with the 



suppression of the number of spin-zero pairs around the 
critical temperature 0, 0] . The important feature is an 
odd-even staggering of the S shape where the S shape of 
the heat capacity for an even number of neutrons is signif- 
icantly morepronounced than the one for an odd number 
of neutrons U • A shell-model analysis [E sug- 
gests that the difference between heat capacities of nu- 
clei with even and odd numbers of neutrons is an impor- 
tant measure for understanding the breaking of neutron 
pairs. Figurejn shows the experimental level densities 
for 93_98 Mo The extrapolated level densities using 
the back-shifted Fermi-gas model shown in Fig. make 
it possible to systematically study pairing correlations 
at finite temperatures in a series of odd- and even-mass 
nuclei. The heat capacities extracted from experimental 
level densities all show a clear S shape as seen in Fig. 
Therefore, the heat capacities of 93 ~ 98 Mo reveal a sys- 
tematic, different from that of other experimental data 
and theoretical calculations 0, El 

The S shape of the heat capacity around the critical 
temperature is not a quantity which reflects details of the 
pairing transition. We cannot distinguish between, e.g., 
quenching of proton-proton (pp) and neutron-neutron 
(nn) pairing correlations from the S shape. In our recent 
papers [9|, llOl lll| , we demonstrated that the suppression 
of pairing correlations around the critical temperature 
appears in thermal odd-even mass differences using the 
three-point formula. For heated systems, the thermal 
odd-even mass difference is the natural extension of the 
odd-even mass difference [13, [13, [lj| observed for nuclear 
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FIG. 1: (Color online) Experimental and calculated level den- 
sities as a function of excitation energy in 93_98 Mo Q. The 
open circles denote the SPA+RPA calculations. The error 
bars show the statistical uncertainties. 
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FIG. 2: (Color online) Extrapolated and calculated level den- 
sities as a function of excitation energy in 93_98 Mo. The open 
circles and triangles denote the SPA+RPA and SPA calcula- 
tions, respectively. 



ground-state masses. 

The aim of this rapid communication is to investigate 
pairing properties in 93_98 Mo. We report on SPA+RPA 
calculations in the monopole pairing model using a de- 
formed Woods-Saxon potential with spin-orbit interac- 
tion and we discuss physical reasons for why the heat 
capacities for Mo isotopes do not show an odd-even stag- 
gering of the S shape. 

Let us start from a monopole pairing Hamiltonian 



H 



(1) 



with 
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where r = n,p and k denotes the time reversed state. 
Here, Sk,r is a single-particle energy and P T is the pair- 
ing operator P T — T Ck, T - Here, we neglect proton- 
neutron pairing interactions because in 93_98 Mo the re- 
spective single-particle orbitals might be too far from 
each other to play a role. By means of the SPA+RPA 
based on the Hubbard-Stratonovich transformation [l( 
the number-parity projected partition function is given 
by 



Z c = Tr \P N Pze 



SPA+RPA 



(2) 



where, instead of the exact number projection, we intro- 
duced the number-parity projections Pjy = (l + ae i7rN )/2 
for neutron number N and Pz = (l + ae™ z )/2 for proton 
number Z (tr denotes the even or odd number parity 
It has recently been shown that in the monopole pairing 
case the SPA+RPA with number-parity projection re- 
produces well exact results 0. The number-parity pro- 
jection is essential to describe thermal properties at low 
temperature. Here, we use the notation u>k,r for the con- 



ventional thermal RPA energies and Afc iT = y c k . 

e' k T = e k .T - Mr - G T /2, and j kiT = e k . T - /x T - X k , T - 
The SPA partition function is obtained by neglecting 
the RPA partition function C^ PA . The thermal en- 
ergy can be calculated from E = —d\nZ c /d(3 where 
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FIG. 3: (Color online) Heat capacities as a function of tem- 
perature T in 93_98 Mo. The solid curve denotes the heat 
capacity extracted from the experimental level density. The 
open squares and open triangles are the calculated heat ca- 
pacities for neutrons and protons, respectively. The total heat 
capacities are represented by the open circles. 



[3 = 1/T. In this work, we use the single-particle en- 
ergies £k, T given by an axially deformed Woods-Saxon 
potential with spin-orbit interaction. We choose the 
Woods-Saxon parameters of where Vq = 51.0 MeV, 
a = 0.67 fm, k = 0.67, A = 20.3, and r = 1.27 fin. We 
have calculated the single-particle spectrum of this po- 
tential, and used it to compute the number-parity pro- 
jected partition function. The deformation parameters 
for the even-even nuclei 94 > 96 ' 98 Mo are estimated from 
the experimental B(E2; 2+ — * 0^ ) values, and are deter- 
mined as /?2 = 0.15, 0.17, and 0.17, respectively. Those 
for the odd nuclei 93 > 95 > 97 Mo are chosen as /?2 = 0.10, 
0.08, and 0.17, respectively. The 25 and 30 doubly de- 
generate single-particle levels with negative energy are 
taken with respect to neutrons and protons outside the 
48 Ca core. The positive energy levels (resonances and 
continuum states) are neglected 0. We adjusted the 
pairing force strengths at G n — 22 /A MeV for neutrons 
and G p — 26.5 /A MeV for protons in order to reproduce 
the three-point odd-even mass differences. 

The level density can be evaluated from an inverse 
Laplace transformation of the partition function Z c in 



TABLE I: Parameters used for the back-shifted Fermi-gas 
level density. 
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the saddle-point approximation as 

Z c e? E 



P(E) 



[2nd 2 InZc/dp 2 } 1 / 2 ' 



(5) 



The calculated level densities as a function of thermal 
energy are shown in Fig. ^ where the calculated level 
density is also plotted for 02 = Q- The calculations re- 
produce very well the experimental data. It should be 
noted that, more precisely, Eq. (JjjJ gives the state density 
and not the level density because the partition function 
Z c in Eq. J2J includes m-degeneracy [l|. This difference 
may have an impact on the canonical heat capacity. 



We now extrapolate the experimental level densities by 
the back-shifted Fermi-gas model of 0, 0] 



Pbsfg{E) = rj 



12V2a 1 / 4 U 5 / 4 a I ' 



(6) 



where the back-shifted energy is U = E—E\ and the spin- 
cutoff parameter 07 is taken as a 2 — 
The level-density parameter a and the parameter E\ are 
given by a = 0.21^ 87 MeV" 1 and E x = d + E peil , 
respectively, where the back-shift is parameterized by 
C\ = — 6.6A~ 032 MeV and the pairing energy impair is 

(3) (3) 

based on pairing-gap parameters A„ and A p evalu- 
ated from odd-even mass differences 01 • The factor 77 is 
introduced in order to reproduce experimental neutron- 
resonance spacings. The parameters used for 93_98 Mo 
are listed in Table [I] Figure [21 displays the extrapolated 
and calculated level densities as a function of excitation 
energy. The level densities from the SPA+RPA calcula- 
tion are in good agreement with the extrapolated ones. 
In Fig. 21 we also compare the SPA+RPA with the pure 
SPA calculation. It is noted, that because of the RPA 
contributions, the calculated level density is significantly 
increased below an excitation energy of ~ 20 MeV. 

To study thermal properties derived from experimental 
level densities, let us start from the partition function in 
the canonical ensemble, i.e., the Laplace transform of the 
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FIG. 4: (Color online) Heat capacities derived using differ- 
ent extrapolations of the experimental level-density curve in 
98 Mo. The extrapolations are performed with a 20% increase 
of (i) the level-density parameter (dashed line) or (ii) the 
pairing-gap parameter (dotted line), respectively. The solid 
curve indicates the heat-capacity curve of Fig. |3] 



1.6 
1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 
1.6 
1.4 

> 1.2 
oj 1.0 

e- 08 

v 0.6 
" 0.4 
< 0.2 
0.0 
1.6 
1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 



" AAA /W<SgS£*£A^ A 93 

a Mo 

A 

A 

A 

□ j-j 

1=1 [ 


AAA AAAfffflffaaYVy^ 

\ 94 Mo 

□ □ □ ni i lllln 2Tj]_^ ^ 

A 

A 

DQ D □ 

a n 

neutron 
a proton 


AAA AA^S&&££t£fA£AA . 95-» 

Mo 

A 

A 

□ □ □ □n3nnnrrx£q- ll - 1| _ | A 

1=1 □ A 

D ° ° □ c 


- AAA A&3SSHKS2AAA. 

* A Mo 

- n n n nr a 


- AAA AA£2SSS!Sf£CAiAA 07, , 

Mo 

A 

A 

□ □□ □n3nrrJ5mtt-i | _ || _ ] A 

□ [ 


Mo 

'' Dn °*\ o 



level density /o(-E'i) 



T (MeV) 

FIG. 5: (Color online) Effective pairing gap as a function of 
temperature. The open squares and open triangles denote the 
effective pairing gaps for neutrons and protons, respectively. 



(7) 



i=0 



where £7, are the excitation energies and 6Ei are the en- 
ergy bins. Then, the thermal energy is expressed as 

oo 

E(Z, N,T) = J2 5Ei P {Ei)Eie- E ^ T /Z(T), (8) 



and the heat capacity is given by 

dE(Z, N, T) 



C{Z,N,T) = 



dT 



(9) 



Formally, the calculations using Eqs. I0)-|0 require an 
infinite summation. However, the experimental level den- 
sities in Fig. gorily cover the excitation energy up to 6.0- 
8.5 MeV. Since the calculated level densities agree well 
with those extrapolated using the back-shifted Fermi-gas 
model of Eq. © , we use this extrapolation to extend the 
experimental level-density data up to an excitation en- 
ergy of —180 MeV before evaluating Eqs. Q-©. The 
5-shape behavior of the heat capacities for 93_98 Mo can 
be reproduced by the calculation, although the calcu- 
lated heat capacities are larger than the experimental 
ones above the critical temperature. The characteristic 
feature is that all of the nuclei 93_98 Mo exhibit a clear 
S shape of the heat capacity around the critical temper- 
ature T c ~ 0.7-0.9 MeV. The calculated heat capacities 



are divided into neutron and proton parts as shown in 
Fig. This figure also indicates that the heat capac- 
ity for the neutron part shows a different behavior from 
the one for the proton part. For the neutron part, the 
critical temperature is — 0.55 MeV and the S shape of 
the heat capacity for an even number of neutrons is more 
pronounced than the one for an odd number of neutrons. 
For the proton part, on the other hand, the heat capaci- 
ties show the same clear S shape around T c ~ 0.75 MeV 
for all of the nuclei 93-98 Mo. Thus, we conclude that 
the experimentally observed S shape of the heat capac- 



ity m 



93-98 



Mo can be attributed mainly to the proton 



contribution. 

We will now discuss how sensitive the extracted heat 
capacities in Fig. |3| are with respect to the level-density 
extrapolation. Figure 01 shows the extracted heat capaci- 
ties using a 20% larger (i) level-density parameter a or (ii) 
pairing-gap parameter -E pa ir- The increased a makes the 
slope of the heat-capacity curve steeper than before; the 
increase of the pairing-gap parameter does not change the 
heat-capacity curve significantly. In either case, the heat 
capacity shows a pronounced S shape which indicates 
that this qualitative feature is very robust with respect to 
changes in the Fermi-gas parameters of the level-density 
extrapolation. 

The S shape has been discussed to be correlated with 
the breaking of nucleon Cooper pairs 0,13- Therefore, 
we further investigate the pairing properties in the calcu- 
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FIG. 6: (Color online) Left: thermal pairing gaps extracted 
from the three-, four-, and five-point indicators of the thermal 
odd-even mass difference as a function of temperature. Right: 
effective pairing gaps for neutrons in the monopole-pairing- 
model calculations of this work. 



lations. Figure [5] shows the effective pairing gap defined 
by 



ld\nZ c 



(3 dG T 



1/2 



(10) 



Around T c ~ 0.6 MeV, A™ ff decreases more rapidly for 
an even number of neutrons than for an odd number of 
neutrons. The suppression of A" ff is well correlated with 
the S shape of the heat capacity for neutrons in Fig. 
13 Furthermore, A^ ff decreases drastically around T c ~ 



0.8 MeV for all of the nuclei 



'Mo. This suppression of 



A^ ff correlates well in temperature with the pronounced 
S shape of the proton heat capacity in Fig. [3] Thus, 
the S shape in Fig. [3| can be understood in terms of the 
suppression of the effective pairing gap. 

The ground-state odd-even mass difference is known 
to be a measure of pairing correlations [l2L fl4| . Extend- 
ing this odd-even mass difference, we have proposed the 
thermal odd-even mass difference as a measure of pair- 
ing correlations at finite temperatures 0, ^| . The three- 
point thermal odd-even mass difference for neutrons is 
given by 



a£\z,n,t)= { -^> 



[B t (Z,N + l,T) 
2B t (Z,N,T)+B t (Z,N-l,T)], (11) 



where the thermal energy B t is defined by B t (Z, N, T) — 
E(Z, N, T) + B(Z, N) and B(N, Z) is the binding en- 
ergy at zero temperature. The thermal odd-even mass 

(31 

difference A n is evaluated using E{Z, N, T) from Eq. 
iJSJ. In our previous work 0] > we extracted the thermal 
odd-even mass difference A^ in 184 W from the experi- 
mental level densities of 183 W, 184 W, and 185 W, and we 

(3) 

obtained a drastic decrease of A„ at the critical temper- 
ature, a signal which can be correlated to a correspond- 
ing S shape of the heat capacity. As we have already 
mentioned in Ref . 0] , the sudden decrease of the ther- 
mal odd-even mass differences is interpreted as a rapid 
breaking of nucleon Cooper pairs. It is therefore inter- 
esting to see whether such a drastic change is found for 
the Mo isotopes as well. 

The upper left panel of Fig. shows the thermal odd- 
even mass differences A^ 3 "* for neutrons as a function of 
temperature. For 94 Mo, A^ starts to decrease rather 
suddenly around T =0.6 MeV, however, this decrease is 
never as drastic as the one observed for 184 W ^3] ■ More- 
over, for 95_97 Mo the thermal odd-even mass differences 

f 3) 

A)j do not display any sudden onset of quenching, rather 
they exhibit a more monotonous decrease with increas- 
ing temperature starting as early as around T =0.2 MeV. 
We now compare these thermal odd-even mass differences 

(3) 

A« to the effective pairing gaps A" ff for neutrons in the 
monopole-pairing-model calculations for 94 ~ 97 Mo (upper 
right panel of Fig. 0). The calculated effective pairing 
gaps decrease suddenly around the critical temperature 
corresponding to the presence of the S shape in the heat- 
capacity curves. At high temperatures of T =1.25 MeV, 
their values are A" ff =0.6-0.7 MeV which is in reasonable 
agreement with the thermal odd-even mass differences for 
94 Mo, however, the thermal odd-even mass differences for 
95 ~ 97 Mo are much smaller than the theoretical estimates 
at these high temperatures. 

It has been demonstrated recently that A^ contains 
additional mean-field contributions when realistic pairing 
forces are used in the calculations , however, the five- 
point odd-even mass difference A n 5 ^ extracts the pairing 
gap. Therefore, we also calculate the four- and five-point 
thermal odd-even mass differences 



A£) (Z n > t ) = \ [M? (Z, N, T) + Ai 3 > (Z,N-1, T) 



(12) 



and 



a£\z > n > t) = ~[aW(z,n + i,t) 
+ 2aw (z,n,t) + a^(z,n- 1,t)1 . (13) 

The lower left panels of Fig.[S]show the thermal odd-even 
mass differences for neutrons a! 4 ' 5 ^ as a function of tem- 



perature, where A„ 



(4) 



and A„ 5) 



treatment as A 

,(4,5) 



are obtained by the same 



The thermal odd-even mass differ- 
ences A^'''' extracted from the experimental data also 



6 



show a monotonous decrease in contrast to the sudden 
quenching of the effective pairing gap around the critical 
temperature in the monopole-pairing-model calculations 
for 94 ~ 97 Mo (lower right panels). 

Thus, the majority of the thermal odd-even mass dif- 
ferences a! 3 ' 4 ' 5 ^ decrease monotonously and do not show 
any sudden quenching, while the theoretical calculations 
suggest the presence of a pairing phase transition at a 
critical temperature of T —0.6 MeV. To possibly explain 
this discrepancy, we would like to point out that the ther- 
mal odd-even mass differences a! 3 ' 4,5 ' are phenomenolog- 
ical indicators extracted from level densities; they may in- 
clude other correlations than just monopole pairing cor- 
relations. E.g., quadrupole correlations, which are not 
taken into account in our calculations, could be impor- 
tant for understanding thermal properties. Such addi- 
tional seniority non-conserving correlations may wash out 
the expected sudden decrease of A„ due to the break- 
ing of nucleon Cooper pairs at the critical temperature. 
Hence, the relation between the odd-even mass differ- 
ences and the pairing correlations at finite temperature 



still remains an open question. 

In conclusion, we have investigated pairing properties 
in the 93_98 Mo isotopic chain of neutron-odd and -even 
nuclei by performing SPA+RPA calculations within a 
monopole pairing model. The calculations reproduce 
very well the experimental level densities, and explain an 
unusual feature recently found for 93_98 Mo, namely that 
the corresponding heat capacities do not show any pro- 
nounced odd-even staggering. The thermal three-, four- 
, and five-point odd-even mass differences in 94-97 Mo, 
which are regarded as measures of pairing correlations at 
finite temperature, were extracted from the experimental 
level densities. They show a monotonous decrease with 
increasing temperature and no sudden drastic quench- 
ing. Further theoretically studies in this direction are in 
progress. On the experimental side, we plan for further 
experiments to also extract proton odd-even mass differ- 
ences at finite temperature. Strong quenching of proton 
pairing correlations as seen on Fig. |S] suggests a dras- 
tic suppression of odd-even mass differences for isotonic 
chains. 
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